OUTP-0619P 



Neutrinos in a Sterile Throat 

Ben GripaioJ*! 

Rudolf Peierls Centre for Theoretical Physics, 
University of Oxford, 1 Keble Rd., Oxford 0X1 3NP, UK and 
Merton College, Oxford 0X1 4JD, UK 
(Dated: December 16, 2009) 

We consider field-theoretic models of a warped extra dimension with multiple 
throats, in which fermions that are singlets of the Standard Model gauge group 
propagate in a separate throat from the Standard Model fields, which we call the 
sterile throat. The singlets mix with Standard Model fields via interactions local- 
ized on the UV brane that connects the two throats. This leads to three, light, 
mostly-active, Majorana neutrinos via a higher-dimensional see-saw mechanism, to- 
gether with Kaluza-Klein towers of mostly-sterile neutrinos, whose scale is set by 
the warp factor in the sterile throat and can be very low if the throat is deep. We 
suggest that a model of this kind may explain all the neutrino data, reconciling the 
LSND result with astrophysical constraints. 
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I. INTRODUCTION 

The observation of neutrino masses and mixing has provided us with the first window onto 
physics beyond the Standard Model (SM). At least inasmuch as the solar and atmospheric 
neutrino experiments are concerned, the data are well described by the SM sirpplemented 
with two or more fermions transforming as singlets under the SM gauge group 1, ^. These 
sterile neutrinos, z/^, interact at the renormalizable level with the three active neutrinos, 
ul, of the Standard Model and the Higgs, H, via a Yukawa interaction (which becomes a 
Dirac mass term for the neutrinos once the Higgs gets a vacuum expectation value) and with 
themselves via a Majorana mass term. Schematically, the SM Lagrangian is supplemented 
with the terms 

6C D XHvlUr + mi'Ri'R + H. c, (1) 

together with the usual kinetic terms for ur- A good fit to the solar and atmospheric 
neutrino data is obtained either by taking the Majorana masses m to be zero (in which 
case the Yukawa couplings A must also be very small and the mass eigenstates are Dirac 
fermions) or by choosing m to be large (around 10^^"^^ GeV, in which case the A are of 
order unity and there are three, light, mostly-active mass eigenstates, with the rest very 
heavy). In the latter, see-saw model, the presence of SM-gauge-singlet fermions at a high 
scale explains why some neutrino masses are so light. 

There is, however, one short-baseline experiment, LSND, that proves to be more prob- 
lematic. The LSND experiment 3|] observes Vf. appearance in a z7^ beam at more than 3(j 
above background, which can be fit by a two-neutrino oscillation with a mass-squared dif- 
ference of around an eV^. If the LSND result and the oscillation hypothesis are correct (the 
MiniBoone experiment, which will independently probe the LSND region, is expected to 
announce its results imminently) then we have evidence for three independent mass-squared 
differences in the low-energy neutrino spectrum, viz. Am^^j ~ 10~^eV^, Am^^^ ~ 10~^eV^ 
and Amj^gj^jj ~ eV^. Naively, these are easily accommodated in the minimal extension of 
the SM with sterile neutrinos and non- vanishing Majorana masses discussed above, by choos- 
ing some of the mostly-sterile mass eigenstates to have mass of order an eV. The oscillations 
relevant for LSND then occur via one or more intermediate, mostly-sterile states. It turns 
out that at least two such states are required for an adequate fit to the data 0, 5, 6, 3,121,0], 
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because of the constraints coming from a number of other short-basehne experiments that 
do not see any evidence for oscillations. 

The objections to this explanation of LSND are both theoretical and observational. From 
a theoretical viewpoint, the explanation of LSND using the minimal extension of the SM 
is unsatisfactory, in that it introduces a light (eV) mass scale in an ad hoc fashion. We 
stress, though, that this objection is purely aesthetic: all values of the fermion masses are 
technically natural, and so assigning arbitrary values to them is no worse than assigning 
arbitrary values to the Yukawa couplings of quarks and charged leptons in the Standard 
Model. Still, an explanation of these values must, ultimately, be forthcoming. 

The observational objection is that extra, light neutrino states are in conflict with stan- 
dard Big Bang cosmology. This objection is more serious, and has become acute following 
the recent measurement of the cosmic microwave background radiation by WMAP [lo| . Ex- 
tra, light neutrino states have two important effects in the early Universe. Firstly, they 
affect the predicted primordial abundance of light elements (in particular ^He) that result 
from Big Bang nucleosynthesis (BBN), because BBN is sensitive to the energy density stored 
in relativistic species, and because BBN depends on reactions involving neutrinos, such as 
z7eP e^n. Using the WMAP data, which provides a precise measurement of the baryon- 
to-photon ratio, the effective number of thermalized, relativistic neutrinos at BBN is given 



by 
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mm 

Nf < 3.4, 95% (two - sided) C. L. 

It thus appears difficult to accommodate even one extra neutrino, let alone two. Secondly, 
neutrinos affect structure formation, since they tend to decouple and free-stream well before 
structure begins to form. The spectrum of initial density perturbations is, therefore, cut- 
off on the smallest scales 15|. The WMAP data, together with observations of large-scale 
structure, then provide a constraint on the energy density stored in neutrinos (more precisely, 
on Vt^h"^) and ergo on the sum of neutrino masses. The constraint is model-dependent; for a 
model with four or five neutrinos, of which three are degenerate, the 95 % confidence level 
bounds are [3] 

Sm^ < 1.4eV, A^"^ = 3 + 1, 

Sm^<2.1eV, A'^ = 4 + l. (2) 
Again, it would appear that multiple neutrino states at an eV or above are only marginally 
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possible. It seems then, that a coherent explanation of all the neutrino data is hard to find.^ 
This state of affairs will become particularly embarrassing if MiniBoone confirms the LSND 
result. 

In this work, we wish to suggest a radical approach to sterile neutrinos. Our starting 
point is the observation that both objections to extra, light sterile neutrino states may be 
surmountable in a model where the extra states are composite degrees of freedom. As far 
as the theoretical objection is concerned, we can explain the existence of an eV scale in 
a natural way by supposing that the required sterile neutrinos are not fundamental fields 
of an effective field theory, valid up to some high energy cut-off, but rather are composite 
bound states of some theory whose couplings grow slowly in the infra-red, and become strong 
around an eV. This is entirely analogous to what happens in the strong interactions, where 
the QCD gauge coupling increases logarithmically in the IR, such that the strong-coupling 
scale is around a GeV and the resulting composite bound states are the hadrons. In this 
way, we find a natural explanation for the existence of a scale that is very light compared 
to other scales in Nature: it arises because the couplings run so slowly.^ 

Moreover, the supposition that sterile neutrinos are composite means that the discussion 
of the effects of the neutrinos on physics in the early Universe is changed. Indeed, at early 
times (in fact, at all energy scales much greater than the compositeness scale), the weakly- 
coupled degrees of freedom of the theory do not consist of the neutrinos per se, but rather the 
degrees of freedom appropriate to the ultra-violet (UV) completion of the theory, whatever 
that may be. 

In the absence of a calculable UV completion, this observation is of limited use. Remark- 
ably, such a completion can be found, by recourse to the celebrated AdS/CFT correspon- 



dence 22|]. In the correspondence, or rather a regularized version thereof 23|, |2J, |25|, the 
UV completion consists of a 4-d conformal field theory of large rank, coupled linearly to dy- 
namical source fields. The theory becomes strongly-coupled in the infra-red (IR), triggering 
spontaneous breaking of the conformal symmetry, and the appearance of a discrete spectrum 
of composite bound states. According to the correspondence, this theory is holographically 
dual to a weakly-coupled (and therefore calculable) 5-d theory formulated on a slice of an 



^ Other approaches to explain the LSND result include CPT- violation \17^ , violation of Lorentz invariance 
18 1, sterile neutrino decay 19 1, and decoherence [iot . 



The idea of compositeness to explain neutrino mass scales was used before in 21 1 
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FIG. 1: The two-throat set-up, with SM fields {H, v,. ■ ■) and sterile neutrinos {^) living in separate 
throats and interacting at a common UV brane. 

AdS^ geometry, that is terminated by branes located at the points z = zq and z = zi in 
the fifth dimension. Going back to the 4-d theory, the extra-dimensional co-ordinate z is 
interpreted as the renormalization group scale, and the branes translate to cutting off the 
CFT in the UV at a scale I/zq, and to strong-coupling occurring in the IR at a scale 1/zi. 
To an observer living in four dimensions, the 5-d fields manifest themselves as Kaluza-Klein 
(KK) towers of 4-d mass eigenstates, the mass spectrum beginning at around 1/zi. These 
KK states are dual to the composite objects that arise when the conformal symmetry is 
spontaneously-broken in the IR. 

Theories of this type have already been extensively invoked in the literature to explain 



the phenomenology of the SM, in particular the naturalness of its electroweak sector 26l. 1271. 
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291 ]. The SM fields are assumed to propagate in a slice of AdS^ (which we call the SM 
throat). By arranging the Higgs field in such a way that its lightest mode is localized towards 
the zi brane, where 1/2:1 ~ TeV, one obtains a solution to the gauge hierarchy problem: 
in the 4-d dual, the Higgs is wholly or predominantly a composite object of the CFT, so 
its mass parameter is insensitive to radiative corrections occurring on distance scales much 
shorter than an inverse TeV. 

Coming back to the issue of sterile neutrinos, we see that if we wish to explain the 
existence of an eV scale in this way, the sterile neutrinos must live in a slice of AdS^ with 
1/zi ~ eV. We immediately see that the SM fields cannot propagate in the bulk of the same 
throat, because there would then be KK excitations of all SM fields beginning at around an 
eV, in clear contradiction with experiment. The sterile neutrinos must therefore live in a 
throat of their own, which we call the sterile throat, and interact with fields in the SM throat 
only through interactions localized on the UV brane that connects the two.'^ This set-up is 



^ The SM fields need not live in a throat at all: they could of course be simply localized on the UV brane 



6 



illustrated in Figure [Tl Such a scenario of field theories in multiple AdS^ throats has been 
recently explored in 30, l3l|, 132], in which multiple slices of AdS^ space are connected via 



a common UV brane. It has a natural dual explanation in terms of separate CFTs that 
become strongly-coupled at their respective IR brane scales and interact with each other 
indirectly via external fields living on the common UV brane. The set-up also provides a 
toy model of recent string theory compactifications 33|| on Calabi-Yau manifolds with fiuxes. 

In what follows, we perform a detailed analysis of the set-up illustrated in Figure [TJ We 
begin, in Section [Tll by developing the theory of a single bulk sterile neutrino propagating 
in a slice of AdS. A gauge-singlet fermion has both Dirac and Majorana mass terms in the 
warped bulk. The general case where neither of these vanishes has received little attention 
in the literature thus far (a brief numerical analysis with sterile neutrinos in the SM throat 
was performed in [m])-*^ There is perhaps good reason for this: we show that the decoupled, 
second-order, differential equations for the extra-dimensional profiles of the sterile neutrino 
KK modes are not of confluent hypergeometric form. They are, however, Fuchsian at the 
origin, and we define the solutions by the two independent series solutions there. Unlike in 
the case of pure Dirac masses 45|, we show that massless or ultralight modes, with tunable 
localization properties, do not occur for natural choices of the boundary conditions on the 
UV and IR branes. Instead, one typically obtains a KK spectrum of Majorana fermions 
with masses beginning at around 

We are most interested in the particular scenario where the bulk sterile neutrino has an 
interaction with the Standard Model fields via a Yukawa-type interaction, localized on the 
UV brane, connecting the active neutrinos, the Higgs and the bulk sterile neutrino.^ We 
calculate, in lll Ct the boundary conditions on the UV brane that follow from this interaction, 
and show that they lead to the presence of an ultralight mass eigenstate in the KK spectrum, 
with mass of order TeV^^o, which is mostly made up of active neutrinos. Plausible values 
for the light neutrino masses result for I/zq in the region of 10^^~^^GeV. The magnitude 
of the ultralight mass is just what one would obtain in the usual see-saw mechanism with 



of the deep throat, with thp_§qji^q_lne;i:a^9h^ pr£|)blem_ap 
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431, 



some other way. 



^ For related references, see 
^ The co-ordinate z will always refer to the sterile throat in what follows. 



We have in mind more sophisticated realizations of the Randall-Sundrum scenario, in which the Higgs i s 
not (5-function localized on the IR brane, but instead propagates in the bulk of the SM throat 
These offer a number of advantages relative to the original Randall-Sundrum example. 
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fundamental sterile neutrinos with mass of 0(1/^0); but results from the mixing of the active 
neutrino with the entire KK tower of sterile neutrinos, and the specific expression for the 
mass reflects this. The spectrum of higher KK states begins at around 1/2:1 and is made up 
of mostly sterile states. At large KK number, the states are separated in mass by "r j z\. 

In Section IIIIl we generalize to more than one bulk sterile neutrino. In this case, the 
analysis is further complicated by the observation that no basis exists, in general, in which 
the bulk Dirac and Majorana mass terms are simultaneously diagonal. This prevents us from 
giving a simple expression for the masses of the ultralight modes. (We do, however, present 
a simple expression for the masses of the ultralight modes in the special case where the mass 
terms can be diagonalized simultaneously.) There are three ultralight modes (one for each 
active neutrino of the SM), predominantly composed of active neutrinos, with order-one 
mixing between the SM flavour eigenstates. As a result of the mixing, neutrino oscillations 
occur between the active neutrino flavour eigenstates. The solar and atmospheric data 
can be explained by oscillations between the mostly-active ultralight modes. We further 
conjecture that the LSND result can be explained by oscillations between the mostly-active 
ultralight states and the mostly-sterile higher KK states. A detailed investigation of this 
requires knowledge of the spectrum of higher KK modes and hence a numerical analysis, 
which is beyond the scope of this work. 

In Section HVt we address the high-energy behaviour of the theory, with the early Universe 
in mind, reverting to the case of a single bulk sterile neutrino for simplicity. The theory is 
well-defined and calculable (via its 5-d gravity dual) at energies up to and beyond the scale 
1/zo (up to the fundamental 5d gravity scale in fact, which we take to be M5 ~ 10^^~^''GeV, 
so as to obtain the correct 4d gravity scale, M| ~ M^/k). 

The crucial consideration, as far as the early Universe is concerned, is whether the extra, 
mostly-sterile states composing the KK tower can come into equilibrium with the primor- 
dial plasma of SM excitations before the epochs of BBN and structure formation. If they 
do (either wholly or partially), then their contribution must be included in the standard 
framework. We envisage two processes by which equilibriation might occur. One is by in- 
teractions in the CFT itself, and the other is by SM interactions. In either case, the actual 
communication between the SM and sterile CFT sectors occurs not at the interaction ver- 
tex, but via the mixing in the propagators, or equivalently in the mass eigenstates. On the 
one hand, we argue that the extra states cannot come into equilibrium by CFT interaction 
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processes. On the 5d side, the explanation for this is that the KK states are locahzed in 
the IR, and can only interact with the SM fields localized on the UV brane through the 5d 
UV-to-IR brane propagator for some field that propagates in the bulk of the sterile throat. 
Now, at large Euclidean four-momenta, these propagators fall off exponentially above the 
KK scale, as e~^^^ On the 4d side, the explanation is that, an Euclidean momenta far above 
l/^i, conformal symmetry is restored, and the composite states are transparent to probes 
on short momentum scales |471] . Thus, even though the number of KK states that are above 
threshold for some process and the available phase space grow as some powers of p, these 
enhancement factors will always be overwhelmed by the exponential suppression at large 
enough p, and interactions between SM fields and the extra neutrino states are negligible at 
energies much above an eV. 

On the other hand, the extra states could equilibriate via SM interactions, just as standard 
4-d sterile neutrinos can 48|]. These processes are not suppressed at high energies. As for 
laboratory oscillations, an understanding of the size of this effect will require numerical 
analysis, which is not carried out here. We hope to perform an analysis of this, and other 
phenomenological effects, in future work. 

Finally, in the Discussion, we address the feasibility of performing a global fit to the 
data and the implications for future oscillation experiments, neutrinoless double-beta decay 
experiments, and astrophysics. 

II. ONE BULK STERILE NEUTRINO 
A. Formalism 

We realize a single throat as a slice of an AdS^ geometry in Poincare co-ordinates^ X'^^ = 
{x'^,z), where z e [zq = l/k,ZQ = e'^^^/k]. The metric (signature mostly-plus) is 



ds^ = j^{r]^,dx^'dx'' + dz^), (3) 



with funfbein = kz6^ . The metric is a solution of General Relativity in a 5-d space with 
negative cosmological constant A5, terminated by 3-branes of positive and negative tension 



^ We use lower-case Greek letters for 4-d spacetime indices, lowercase Latin for 3-d space indices and 
uppercase Latin for 5-d spacetime indices. 
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located at Zo and zi, respectively. The AdS curvature scale is given by 



2 ^ U) 



where M5 is the 5-d gravity scale. The brane tensions, A4, are given by 

k=^,- (5) 
6M| ^ ' 

For a multi-throat geometry jsi] , we splice together slices of AdS (with fifth-dimensional 
co-ordinates z,y,...) by identifying the point zq = 1/k with yo = l/l and so on, corre- 
sponding to a common UV brane. The curvature scales k,l, . . . are set by the values that 
the cosmological constant takes in the bulk of each throat, as per Applying the Israel 
junction condition on the UV brane leads to the relation 

* + ' + — i^. (6) 

The five- dimensional Dirac algebra, {7"^, 7^} = 21]^^, is represented by 



where = (l,cr'), = (1,— cr*), and a* are the usual Pauli matrices. The 5-d charge 
conjugation matrix 

satisfies the relation 

C57^C5-^ = +(7^)*, (9) 

where t denotes the transpose. 





6M|' 













"J' 












B. The Bulk Action and Equations of Motion 



The smallest irreducible representation of the 5-d Dirac algebra is carried by a four- 
component Dirac spinor \E'. The most general Hermitian and 5-d Lorentz-invariant action 
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for is then 



S^- j d^xdz (^^r*^DM* - cA;** - (iA;**C5^) + H. c, (10) 

where = e^'-y^ are curved-space gamma matrices and Dm = du + includes the spin 
connection. The parameters c and d are of order unity in a natural theory and we may take 
c to be real, without loss of generality. In terms of the re-scaled Weyl spinors, -^aand x", 
defined such that 

* = {kzf {A , (11) 

the action is 

+ ^(xV' + V'X) + ^(^^ - xx) + ^(^V' - xx)) ■ (12) 

The bulk equations of motion are 



_ _ c_ d— 

= -ia ■ dip - d^x + -X + -V', 

z z 

c d 

^ -ia ■ dx - d,ijj + -ijj - -x- (13) 

z z 

In general, the 4-d mass eigenstates will correspond to Majorana fermions, Cn{x), satisfying 

-ia ■ d^n + rriniri = 0, (14) 
such that the Kaluza-Klein expansion takes the form 

X = ^n9n{z)in{x). (15) 

The bulk equations of motion then become 

= Zg'^ - CQn - dfn + mnzfn. (16) 

In the case of pure-Dirac bulk masses {d — 0), these equations are easily decoupled; jnl \fz 
and gnl \fz are given by linear combinations of the two independent Bessel functions of order 
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c + 1/2 and c — 1/2, respectively. The case c? 7^ is not so straightforward, as we illustrate 
by considering the special case where we take d to be real. Then, differentiating the second 
equation of (|T6l) and substituting for Qn and g' we obtain the following uncoupled equation 
for /n, 

^^^\z z-d/mj^" V + z{z-d/m^) ) ' ^ ^ 

with a similar equation for 5f„. This has regular singular points at z = and z = d/rrin in 
z eC, and an irregular singular point at 2; = cxd (as one can show by substituting z = 1/^). 
Thus, (fT7|) is not Fuchsian, nor is it of confluent hypergeometric form, except when c? — >■ 
(in which case the two regular singular points coalesce in a single regular singular point). 

Given that the origin is a regular singular point, we are able to find two independent series 
solutions about this point by the method of Frobenius. These serve as suitable definitions 
for the independent solutions, at least for 1^1 < d/rrin. This will also turn out to be the 
most useful definition for phenomenological reasons, since the leading order behaviour of the 
solutions near the orgin determines the masses of ultralight modes, with rrinZi <^ mnZo ^ 1. 
Returning to the most general case where d e C, we make the ansdtze 

f — V /,{")^"+™ ^ _ V hW ■y(^+'^ ('\'A\ 

Jn — ^?Ti=0"nm^ ' i/™ ~ ^m=0'^nm'^ ' \^°) 

that yield 

= S„=o ((« + m + 0)4"^^"+'" + dhZz'^-'"^ - m„6(« ^^+"^+1) , 

= S„=o ((/5 + m - c)6(fj,z'^+'" + rfai"^;2"+'" + m„a(-^z°+™+i) , (19) 

upon substitution. The indicial equation, at m = 0, has a solution with non-vanishing 
and 6^0^ iff. (3 = a, and 

a' = cc + dd ^ 6i? = -^ai"o^ (20) 
The two roots a = a± are, real, equal in magnitude and opposite in sign. Let us re-define 



a 



+ Vcc + dd = a; (21) 



Then, the two roots are given by a± = ±a. The subsequent unknown coefficients ai^>o 
and &im>o determined iteratively by the analogue of (fT9l) . but containing a linear super- 
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position of the series for both roots a±, viz. 

= S^=oS„=«± ((a + m + c)a("iz-+™ + - m„6(r^^^+™+i) , 

= S^=oSa=a± + m - c)6(f2'S''+'" + ^il^™ + m„a(fiz"+'"+i) . (22) 

The coefficients of the two series with differing roots a± are determined independently of 
each other. Thus we obtain two independent solutions for all c and d (with the exception 
of the special case where \/ + dd is half-integral, and the second independent solution 
contains a logarithm). We can use these series to define the independent solutions as far as 
the singularity at \z\ = |(i/m„|, and presumably by analytic continuation beyond. 

We can also solve the bulk equations of motion (fT6|) at large values of \mn\z. At leading 
order, they become 

g'n ^nfn, (23) 

with solution 

fn ~ Az^ sin |m„|z + Bz^ cos \mn\z, 

Qn ~ —Bz^ sin \mn\z + Az'' cos \mn\z. (24) 

At next-to-leading order, the equations (fT6|) have consistent solution only if p = 0, whence 
the asymptotic solutions at large argument are 

fn^A sin \mn\z + B cos 

^ —B sin \mn\z + Acqs \mn\z. (25) 

These can be matched onto the solutions at small argument, given a suitable analytic con- 
tinuation. 

C. Boundary Dynamics 

We wish to consider the boundary conditions a.t z = zq and z = Zi that result from 
dynamics localized on the UV and IR branes. In particular, we shall be interested in the 
boundary conditions at 2; = 2;o that result from the UV-brane-localized interactions between 
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the fermions in the sterile throat and fields in the SM throat. This is most conveniently done 



following the formalism of 



491 ] ■ To wit, we consider the theory compactified on the interval 



z e [zq, zi], as opposed to an orbifold. The bulk equations of motion ( fT3l) are obtained from 
the vanishing of the variation of the action ( |T2i) . which also includes the boundary variation 



5S 



(26) 



This boundary variation vanishes iff. 



ai z = zq, zi, where 



(27) 



(28) 



(M^ and A'^ may include boundary derivatives, and the ± signs account for partial inte- 
grations of these.) What are the particular forms of the general boundary condition fl27|) 
that result from adding brane-localized terms to the action? At least in simple cases, this 



question can be answered unambiguously 49|] by first adding the localized terms in the bulk 
at a distance e away from the relevant brane, viz. dX z = Zq = z^ + t oi aX z = z^ = Zi — t and 
solving the jump equations there. Then, one matches the solution to the original boundary 
conditions aX z = zq or zi and takes the limit e ^ 0. This slightly convoluted procedure 
avoids dealing with functions that are discontinuous at the boundary. 

As an example, let us consider how the standard boundary condition xl^^ = is modified 
by a Dirac mass term mixing the field ip and a brane-localized fermion rj{x).^ We add to 
the action f[T^ the UV-localized term 



5S 



(29) 



which is manifestly Hermitian and 4-d Lorentz- invariant. Later on, we will take t] to be 



In the holographic dual 45|, tp{x, zq) represents a dynamical source field, external to the CFT, that mixes 
with the CFT and with the external field ri{x). 
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some linear combination of the SM neutrinos. The new bulk equations of motion are 








d- 



—i(7 ■ dip — dzX H — X H — ^ + ^(^ ~ ZQ )hk^^'^r] 



z 



-ia ■ dx - d^ip + -i>- -X, 
z z 



= 6{z - z+){-ia ■ dr] + bk^/^ip). 
Assuming that rj and are continuous at Zq , the jump equations are 

= -ta.d[xU + [dzi:U, 



(30) 



(31) 



where [x]z+ denotes the jump in % at z = Zq . Matching to the original boundary condition 
x\zo = 0, and taking the limit e — 0, these reduce to the second bulk equation of motion 
([T3|) evaluated at the boundary, together with the new boundary condition 



= {-X + bk^^^v) 



(32) 



We may use the third equation of ( l30l) to write this condition purely in terms of bulk fields 

as 



= [—ia ■ dx + bbkil)) 
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(33) 



Since the bulk fermions are gauge singlets, we are also free to add Majorana mass terms, 
localized on either the UV or IR branes, of the form 



6S = — J dt^x {atpip + H. c. 
which leads to a modified boundary condition of the form 

= (-X + ai/j) 



2(^,2 J 



(34) 



(35) 



Such terms present no added difficulty in terms of computation, but nor do they change the 
character of the results. We disregard them in what follows for the sake of clarity. 



D. The Ultralight Mode 



In the case of pure-Dirac bulk masses, d = 0, the boundary conditions x\ 



for 



^1 



Z=20,2l 



0) result in a massless zero mode. No such massless mode is obtained with d ^ 0, 



15 



as we now show. For mo = 0, the general solution to the bulk equations of motion (fT6l) is 
given exactly by the leading order part of the series solution f|T8|) . which we write as 



/o = a^z"' + a_z 

a + c_ ^ -a + c_ _^ 

go = — ^^+^ ^ — "-^ • ^^^) 

A non-trivial solution with one of a± non-vanishing, and satisfying the boundary condition 
x\z=zo,zi =^ go\z=zo,zi = is possible only ii a = ±c ^ d = 0, implying that zero modes 
are only obtained in the pure-Dirac limit. The same argument shows that no modes with 
masses rrinZi ^ 1 are present given such boundary conditions, and we conclude that the 
spectrum begins with modes whose mass is at least of order 1/zi. 

What happens if we add a UV-brane-localized fermion 77? The boundary condition at 
z = zq is changed to (1551) . or equivalently 

= -Wngnizo) + bbkfn{zo), (37) 

together with the original boundary condition at z = zi, namely 

= gn{zi). (38) 

Now, as 6 0, (137|) is trivially solved by mo = 0, so that we are guaranteed the existence 
of a massless mode. This is no surprise: as 6 ^ 0, the fermion t] becomes decoupled from 
the bulk fermion ijj and is manifestly massless. Similarly, for sufficiently small values of b, 
we expect to find an ultralight mode in the spectrum with mo-Zi <^ 1. To find the mass of 
this mode, we can use the small argument solutions of fl36|l . in terms of which the modified 
boundary conditions become 

= -^((a + c)a+ZQ + {—a + c)a^ZQ°') + bbk{a+z^ + a_ZQ"), 

(Jj 

„ a + c_ ^ -a + c_ _^ , . 

= — a+z" — a-Zi . (39) 

After some algebra, we find that these equations have non-trivial solutions iff. mo is such 
that 

\mo\^ = (40) 
c — a 

This expression is manifestly real and positive. We find in either case that the mass is of 
order |6p/c, so is ultralight only if 6^ <^ zq/zi. We shall learn in the sequel that b naturally 
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takes values of order TeV/fc, so we get an ultralight mode if z^^ ^ (TeV)^/fc. Note also 
that the relative fractions of the bulk and brane-localized fields that make up the ultralight 
mode are -j/^ixir^^foifoil, so that the ultralight mode is mostly active and the higher 
KK modes will be mostly sterile. 

Can we say anything more about the spectrum of higher KK modes? For low values of n, 
we would need to solve the bulk equations of motion ( JT6l) with given boundary conditions 
numerically to find the spectrum. However, at larger KK numbers, such that 2;f ^ -C m„ <^ 
Zq'^, we can approximate the solutions to (JT6l) in the region of the IR brane by the large 
argument solutions (l25|l . namely 

/„ ~ A sin |m„|z + B cos 

~ — i? sin |m„|2; + A cos Im^l^. (41) 

Now, a general boundary condition on the IR brane of the form {x + CLip)\zi = 0, where a is 
a constant, will lead to an expression for the ratio of the unknowns A/B as a rational linear 
function of tan |m„|2;i. The allowed values of |m„| are those for which this ratio coincides 
with the ratio of the unknowns on the UV brane, which will be given by some function 
of |m„|2;o, /(|m„|2;o) say. Imagine that |m„| is such a value. If |m„| is increased by ir/zi, 
the rational function returns to the same value, whereas /(|m„|2;o) changes by a negligible 
amount: /(|m„|2;o) /(|m„|zo) + ^f'{\fnn\zo)- Thus, up to corrections of 0{zq/zi), we 
see that KK modes at large n are separated in mass by an amount tt/zi. 

III. MULTIPLE BULK STERILE NEUTRINOS 

With just one bulk sterile neutrino, we obtain only one massive, but ultralight, mode, 
which is mostly composed of a linear combination of the active neutrinos. However, the 
solar and atmospheric neutrino data indicate two independent mass-squared differences, so a 
realistic sterile throat model should contain at least two bulk fermions in the sterile throat.^ 
With this in mind, we now generalize the formalism to the case of multiple bulk sterile 



^ Or, alternatively, more than one sterile throat. 
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neutrinos coupled to multiple brane- localized fermions 77^, the bulk action generalizes to 

..-/,..(-,>.^-.,...,.l(v.,,«-,„-.v,) 

where the matrices Cij and dij are Hermitian and symmetric, respectively. To this we add 
the boundary action 

(55 = - y d^x (^-iri^a ■ df]^ + bjik^/^rj^tP' + blj^rj^^ ^. (42) 

where bji is a complex matrix. The Hermitian matrix Cij can be diagonalized by a unitary 
transformation that acts identically on ip^ and but the symmetric matrix d^j will not 
in general be diagonal in this basis. Similarly, if we choose a basis in which the masses of 
charged leptons are diagonal, then we must also consider a general form for the matrix bji. 
As in the case of one bulk singlet fermion, we solve for the mass eigenstates by expanding 

# = E„/;(^)a(x), 

= S„^^(z)e„(x). (43) 



A series solution for the bulk equations of motion generalizing ([T6|) is then given by 

fi — y ^^Wi^QA+m 

gl, = S^,,6W^z^+-, (44) 



where a\ are the eigenvalues of the matrix 



d'^ 



(45) 



The boundary conditions are, at z = Zi, 

Xi.. =0, (46) 



and, at 2 = ^0? 



0. (47) 
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Substituting from (jH]), the determination of the masses of the ultrahght modes is simply 
an exercise in hnear algebra, albeit somewhat involved in the general case. To illustrate the 
character of the result, we find the masses in the simpler case where the mass matrices c 
and d are simultaneously diagonalizable. Then the bulk equations of motion for different i 
decouple and the solutions at small arguments, satisfying the IR boundary condition, are 
approximately 



The UV boundary condition then reads 



c* + a* zi c* — a* zi 
9n = ^{ i-r - • (48) 



where 



m„a* = Bija\ (49) 



Now, since 



m„a* = Bija^ rrina^ = Bijo' (51) 

we find that |m„p is an eigenvalue of the matrix BB. By considering the 6 — limit, we see 
that there will be one ultralight mode for each brane fermion, each with mass of 0{b'^/zo), 
where 6 is a typical element of the matrix bn. The ultralight modes are predominantly 
composed of active neutrinos, with 0(1) mixing between the flavour eigenstates. 



A. Oscillation Phenomena 



The pattern of neutrino masses and mixings that arises for two or more bulk sterile 
neutrinos is qualitatively similar to that which is observed in so-called 3 + 2 models. These 
have three ultralight modes, which are mostly-active, with 0(1) mixings between the three 
active neutrino flavours of the SM. On top of this, there are two additional, light neutrino 
states, which are mostly sterile. 

In the 3 + 2 models, it is argued that this structure gives a good fit to the oscillation data, 
with solar and atmospheric data explained by oscillations between the ultralight modes, and 
the LSND result explained by oscillations occurring via the mostly-sterile states. 
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The sterile throat models are similar, with the major difference that there are now a 
huge number of mostly-sterile states (the higher KK modes). The beam energy at LSND is 
roughly an MeV, meaning that around 10® mass eigenstates in the KK tower are available for 
oscillation. Because the KK masses are not quite commensurate, the oscillation signal will 

n 

appear to be essentially aperiodic [36'], but this is of little relevance for LSND, which does 
not resolve the oscillations. Indeed, roughly-speaking, the only criterion that LSND imposes 
on the mass scale (or scales) is that it be large enough such that significant oscillations occur 
over the short baseline. This is guaranteed by choosing a large-enough KK scale. 

Can a sterile throat model explain the spectrum of solar and atmospheric mass-squared 
differences? We found above that the ultralight modes have masses of order h"^ / zq. The 
coefficient h comes from the Dirac mass term localized on the UV brane and its size, assuming 
Yukawa couplings of order unity, is set by the expectation value of the Higgs field on the 
UV brane. 

Now, in the most realistic models of warped electroweak symmetry breaking 0, , the 
4d Higgs field corresponds to the zero mode of the fifth component of a bulk gauge field. 
This set-up is holographically dual to realizing the Higgs as a pseudo-Goldstone boson in 4d, 
and protects the Higgs mass parameter from dangerous radiative corrections down to the 
KK scale. The zero mode has an IR-localized profile A^{y^x) ~ yA^{x), where y is the fifth 
co-ordinate in the SM throat. At the UV brane, the Higgs vev is therefore expontentially 
suppressed relative to its vev on the SM IR brane (which is k in 5d units) and h ~TeV.^° 

Alternatively, if the SM is simply taken to be localized on the UV brane, we see directly 
that h has to be of order TeV, so as to get the correct masses for SM fields. This value is, 
of course, unnatural though. 

Thus to get the right masses for the ultralight modes, we need to choose k = I/zq ~ 
2Qi4-i5QgY Normally, when one considers warped compactifications, one simply assumes 
that all fundamental 5d mass scales are roughly the same. If we were to apply this logic here, 
we would conclude that the 4d Planck scale, M4 would be too low. However, we argue that 
not only is a small hierarchy between k and M4 desirable here, but also that it is essential for 
the consistency of the 5d effective field theory formalism that we have employed throughout. 
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In models where the Higgs is the hght mode of a bulk scalar 46[ , the extra-dimensional profile of the bulk 



Higgs, and ergo the vev on the UV brane, can be tuned. This gives further freedom in the value of b. 
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Indeed, the 5d effective field tfieory is valid up to the scale at which 5d gravity becomes 
strongly coupled, M5. In order to be able to trust the warped geometry solutions of the 5d 
Einstein field equations, we require that the bulk cosmological constant and brane tensions, 
Ag^^ and iA^"', should be much less than M5. Since the curvature scale k and the 4d Planck 
scale are given by 



respectively, we see that a hierarchy between the two is obhgatory, and that A5' /M5 ~ 
10"^^^ or so suffices to reproduce the correct Planck scale. Wc thus conclude that masses 
for the ultralight modes in the range consistent with solar and atmospheric oscillation data 
can be easily achieved. These oscillations then occur via the order-one mixing between the 
ultrahght modes. 



We have seen above that sterile neutrinos in a separate throat to the SM may provide an 
explanation for all neutrino oscillation data, with natural explanations for light and ultralight 
masses in terms of a high scale, via the mechanisms of warping and a higher-dimensional 
see-saw. We now consider the effect of the extra KK towers of mostly-sterile neutrinos on 
physics at high energies, in particular in the early Universe. 

As far as CFT interactions are concerned, a key observation is that, as we show below, the 
extra states (the higher KK modes) are all localized in the IR. This means that interactions 
with SM fields (which occur only on the UV brane) can proceed only via the 5d UV-to-lR- 
brane propagator for some bulk field that interacts with the mostly-sterile higher KK fermion 
modes. Examples of such fields include the bulk sterile neutrinos themselves, the graviton, or 

Note that the cut-ofFs of the 5d theory (M5) and its 4d dual (fc) need not be and are not the same. The 
correspondence is only valid up k, because at shorter scales the extra dimensions of the AdS^ x ^5 (or 
similar) manifold show up, and the simple 5d/4d correspondence breaks down. Since the strongly-coupled 
CFT is defined by its weakly-coupled dual, the cut-off of the 4d theory is necessarily k. In the multi-throat 
case, the correspondence breaks down at the smallest curvature scale. 




(52) 



IV. HIGH-ENERGY BEHAVIOUR 
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some other SM-gauge- singlet field that propagates in the sterile throat. A generic property 
of such propagators in a slice of AdS is that they decay exponentially for values of the 
Euclidean four- momentum far above the IR scale 1/2:1. We will show this explicitly for bulk 
sterile neutrinos with non-zero bulk Majorana masses, subject to the boundary conditions 
considered previously. It has already been shown for other bulk fields in js^. This means 
that equilibriation processes involving interactions in the CFT are negligible at energies 
and momenta much above an eV. In the 4-d dual, the explanation for this fall-off is that 
such interactions proceed through the strong dynamics of the CFT, which involves a form 
factor. The CFT in transparent to probes on scales much shorter than an inverse-eV, which 
accounts for the fall-off of the form factor. 

This conclusion does not hold for SM interaction processes, in which equilibrium is at- 
tained via pure SM interactions and transmitted to the sterile states via oscillations. These 
oscillations occur via the two-point 'interactions' corresponding to the couplings between 
the SM neutrinos and the elementary source field (dual to ip{zo, x)), and between the source 
and the CFT. These two-point couplings do not go through the strong dynamics of the CFT 
and do not come with a form factor. 



A. Localization Properties of the Modes 

To show that the higher KK modes are localized in the IR, we show that they are 
normalizable even as the UV brane is removed, that is, in the limit as 2:0 ~^ 0. We consider 
the case of a single bulk sterile neutrino for simplicity. In order to canonically normalize the 
kinetic terms for ^n{x), we see from ( IT2|) . (ITSl) and (!29l) that 

1 = ±.\g„(zo)\' + Hdz {\Uz)\' + \g4z)\'), (53) 

bOK J zo 

where fn{z) and gn{z) are solutions of (fT6l) satisfying the boundary conditions (1371) and 
fl'n(^i) = 0. For modes such that m„ remains finite a.s Zq —>■ (namely all modes with finite 
n bar the ultralight mode with n = 0), the boundary conditions reduce to gn{zi) = and 
fn{zo) = 0, and the normalization integral ( l53l) reduces to 

"dz (|/„(z)r + K(z)r). (54) 
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Near z = Zq^wq can approximate the solutions /„ and Qn by their values at small arguments. 
They thus take the form 

^(0"+0i• 
\ d zq d zq J 
Substituting these solutions into fl54l) . it is clear that the integral is finite as zq 0. There- 
fore all modes with n 7^ are normalizable as 2:0 ^ and are localized in the JR. We repeat 
that this conclusion does not hold for the ultralight mode with n = 0, for which uiqZq is 
finite as 20 ^ 0.-^^ 



B. The Brane-to-Brane Propagator 

We now calculate the UV-to-IR brane 5d-propagator for the bulk sterile neutrino at large 
4-momenta, showing that it is exponentially suppressed for Euclidean momenta above the IR 
scale 1/ Zi. At large 4-momenta, the bulk equations of motion (|T3l) can be self-consistently 
truncated as 

Q = ia ■ dip + dzXi 

Q = -ia-x + d,ilJ. (56) 
We therefore define the Green's functions at large Minkowskian 4-momenta as 

Sp, (57) 



subject to the boundary conditions 

S^{z,,z',p) = S^^{zuz',p) = 



—Sayizo,z',p) = a^.p^,S^^{zQ,z,p) + — 

Zq ^ Zq 

The Green's functions defined in (1571) are given by 



cr'^0Pt,S^{zo,z',p) + —Sa^{zo,z',p) = a^^p^S'^^{zo,z',p) + ^S2{zo,z',p) = 0. (58) 



(si 




H 









(59) 



^"^ The ultralight mode is, in general, localizable nowhere. This is similar to the zero mode of a bulk gauge 
field, or a Dirac fermion with c = 1/2. 
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where Gi^2 are themselves Green's function solutions of 

{dl-p')G^,2{z,z',p)^5{z-z'). (60) 

After some tedious algebra (which we spare the reader), we obtain the UV-to-IR brane 
propagators for spacelike Minkowskian momenta 
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where C = Performing the analytic continuation po ^ ^Po? we see that the propa- 

gators do indeed fall off exponentially for Euclidean 4-momenta greater that the IR scale 



C. High-Energy Phenomenology 

The exponential fall-off of the bulk propagators essentially means that CFT interaction 
processes involving the higher KK states and SM fields have negligible cross-section at 
energies much above an eV. Indeed, one can imagine that there arc enhancement factors, 
such as mass threshold and phase space effects, which grow as powers of the energy, but 
these will always be overwhelmed by the exponential suppression of the propagator at large 
energy. 

In particular, we see that the extra states can never come into thermal equilibrium with 
SM fields in the early Universe through processes of this kind, at least when the temperature 
of SM fields is above an eV. By the time the temperature reaches an cV and the extra 
states have a chance of reaching equilibrium, nucleosynthesis will be complete, and structure 
formation will already have begun. 

The extra states could however, come into equilibrium via SM interaction processes. In 
order to assess to what extent the extra KK states reach equilibrium in this way, one would 
need to perform a numerical solution of the Boltzmann equations in the early Universe, 
taking into account the effects of neutrino mixing and the multiple neutrino mass thresholds. 

We make just one remark in this regard, which is that even if there are multiple massive 
states in complete equilibrium with SM fields, they are not necessarily in contradiction with 
WMAP data. This is because increases in either YuyVriy or can be compensated in the 
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data fit, by moving the present Hubble parameter in opposite directions. Thus, increasing 
both simultaneously has a compensatory effect, reflected for example in ([2]) [l6|. 

V. DISCUSSION 

In the foregoing, we have examined the impact of SM-gauge-singlet fermions (a.k.a. sterile 
neutrinos) propagating in the bulk of a separate throat on physics in the SM throat. We 
find that the bulk sterile neutrinos give rise to three, mostly-active, ultralight Majorana 
neutrinos in 4-d, with mass of order TeV^2;o; where 1/2:0 is the scale of the UV brane, or 
equivalently the cut-off of the dual CFT. The bulk sterile neutrinos also lead to KK towers 
of mostly-sterile Majorana modes, beginning at the IR brane scale, or equivalently the IR 
strong-coupling scale of the CFT. There is order-one mixing in the mass eigenstates between 
the SM flavour eigenstates. What is more, these conclusions hold for all order-one values of 
the dimensionless mass parameters in the bulk or on the boundary: the phenomenology is 
quite generic. 

We have suggested that such a set-up could plausibly explain the entirety of the existing 
neutrino data, including LSND and cosmological observations. We remark in passing that 
there is another experiment, the Heidelberg-Moscow experiment, which in one analysis 51 1 
claims a signal for neutrinoless double-beta decay, interpreted as evidence for Majorana 
neutrino masses. The measured effective neutrino mass |Sjymjy?7g^|, where Ueu is the mixing 
matrix element between the electron-neutrino and the mass eigenstates, is given by the 



ranges 0.11 - 0.56 eV [51i or 0.4 - 1.3 eV 



52l |. depending on the values one takes for the 



relevant nuclear matrix elements. These values could also be explained by a sterile throat 
scenario, with multiple states of eV mass and small matrix elements. 

The next step in the analysis of sterile-throat models is a detailed numerical fit to the 
data. Is this feasible? As regards LSND, or physics in the early Universe at MeV energies, 
this in principle involves finding the masses and mixing matrix elements of around 10^ KK 
modes, given a spectrum beginning at around an eV. However, the active component of the 
KK mass eigenstates goes to zero at large KK number, so we expect that including only a 
relatively small number of modes will prove to be a good approximation in practice. 

Supposing that a sterile throat model does fit the existing data, how might we find 
conclusive evidence for it? The oustanding signature of an extra-dimensional explanation 
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of LSND is the KK tower beginning, at low energy. This has a dramatic effect on the 



36| . Whether any neutrino 



periodicity of oscillations, which appear essentially aperiodic 
beam experiment will ever be able to resolve the oscillations sufficiently to observe this 
seems rather doubtful, but the observation of even a few oscillation components with mass- 
squared differences of similar order would already be a strong hint. Even if the presence of 
a KK tower in the neutrino oscillation data could be inferred, we would still be faced with 
the question of the nature of the relevant extra dimension: would it be warped, or flat, for 
example? These two scenarios are easily distinguished if the the masses of the flrst few KK 
states can be measured with reasonable precision, as they show integer spacing in the flat 



case 
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36|, but do not in the warped case. 
Another possible effect of the KK tower is on high-energy astrophysical processes, for 
example invisible energy loss in stars and supernovae. Such objects have temperatures cor- 
responding to a large number of KK states being above threshold, albeit with tiny individual 
couplings. It would be very interesting to explore the collective effects on such processes of 
a sterile throat model. Again though, detailed numerical simulations are required. We hope 
to address this in future work. 

Finally, we remark that the idea of an additional sterile throat, but with other bulk flelds. 



may lead to other interesting phenomenology. As examples, we suggest the axion 53(] or 
quintessence flelds. 
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